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Abstract 


The transverse instability of a bunch in a circular accelerator is elaborated in this paper. A new 
tree-modes model is proposed and developed to describe the most unstable modes of the bunch. 
This simple and flexible model includes chromaticity and space charge, and can be used with any 
bunch and wake forms. The dispersion equation for the bunch eigentunes is obtained in form of 
a third-order algebraic equation. The known head-tail and TMCI modes appear as the limiting 
cases which are distinctly bounded at zero chromaticity only. It is shown that the instability 
parameters depend only slightly on the bunch model but they are rather sensitive to the wake 
shape. In particular, space charge effects are investigated in the paper and it is shown that their 
influence depends on sign of wake field enhancing the bunch stability if the wake is negative. The 
resistive wall wake is considered in detail including a comparison of single and collective effects. A 
comparison of the results with earlier publications is carried out. 
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INTRODUCTION 


Two kinds of single bunch transverse instability in circular accelerators are distinguished 
at present: the head-tail one and the transverse mode coupling instability (TMCI). First of 
them has been studied hrst by Pellegrini |l| and Sands who were treating the transverse 
bunch oscillations as a combination of uncoupled multipoles oc exp{im(j)) with 0 as syn¬ 
chrotron phase. Role of synchrotron amplitude has been investigated later by Sacherer who 


was introducing the conception of radial modes A 


ot of papers dealing with this were 
8|]). The most important conclusion 


published later, including reviews and books (see e.g. [?, 
is that the head-tail instability is possible at non-zero chromaticity only and depends on its 
sign. 

Other type of instability was observed in many electron machines regardless of chro¬ 
maticity and being primary referred as “transverse turbulence”. Its hrst explanation has 
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been proposed by Kohaupt who used a simple model of two particles propelling each other 
through the mediation of constant wake helds jd]. Later the theory has been evolved on 
the base of Vlasov equation treating the instability as a result of reciprocal influence of 
neighboring pairs of the head-tail modes which frequencies approach each other due to a 


wake held, whence the term TMCI appears [5H8[. Chromaticity has been ignored in these 
investigations as a factor of secondary importance. Some space charge effects has been con¬ 
sidered later for both head-tail and TMCI instabilities jo-Ul. However, specihc models like 


the hollow bunch in the square potential well have been used in {q], and only high space 
charge limit was investigated in Ref. Blul. 

An important observation is that, in all considered cases, the lowest multipoles m = 0, ±1 
or their combinations are the most prone to instability and pose a major threat to the 
accelerator operation. It is shown in this paper that the general theory can be essentially 
simplihed and unihed being restricted by these cases. 

General 3-modes model is developed for this, providing simple analytical formulae for 
all types of single bunch instability. The model is built on the base of exact solutions 
for hollow bunch which are obtained in the paper without any additional assumptions like 
truncated multipole expansion. The dangerous modes are separated and extended to realistic 
distributions with arbitrary wake, chromaticity and space charge included. 

It is shown that the instability parameters depend only slightly on the bunch shape being 
rather sensitive to the wake function form. Resistive wall wake is investigated especially to 
assign which effect is more dangerous in specihc cases: either TMCI or collective mode 
instability. Chromaticity is taken into account in all the cases demonstrating an absence of 
clearly dehned boundary between the head-tail instability (separated multipoles) and TMCI 
(coalesced multipoles). Space charge effects are studied with arbitrary ratio of betatron tune 
shift to synchrotron tune. In contrast with earlier papers Is il], the investigation is not 
restricted by simplihed bunch models or by limiting cases only. It is shown that the space 
charge ehect can be either advantageous or not for the bunch stability, dependent on the 
wake sign. 
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II. BASIC EQUATIONS AND DEFINITIONS [111] 


Linear synchrotron oscillations are considered in this paper being characterized by am¬ 
plitude A and phase 0, or by corresponding Cartesian coordinates: 


0 = A cos 0, 


M = A sin 0. 


( 1 ) 


Thus 0 is azimuthal deviation of the particle from the bunch center in the rest frame whereas 
the variable u is proportional to the momentum deviation with respect to the central bunch 
momentum (proportionality coefficient is not a factor in this paper). Steady state of a bunch 
will be described by the distribution function F{A) and by corresponding linear density 


p{0) = / F{V0^ + u^) du 

J —OO 

with the normalization conditions 

f*oo /*oo 

2n F{A)AdA = l, / p{0)d0 = l. 

Jo J—CO 


( 2 ) 


(3) 


It is convenient to present coherent transverse displacement of the bunch in the point of 
longitudinal phase space {A, 0) in the form [ll| 


Y (A, 0) exp [ - i{Qo + () 0 - i {Qo + ly) ^ot ] 


(4) 


where Qq is the revolution frequency, Qq is the central betatron tune, u is an addition due 
to a wake held, and ( is the normalized chromaticity: 

^oQp _ ^ 


c = 


Q' 

p 


(5) 


1 / 7 ^ — a ’ 

with ^ as usual chromaticity, and a as the momentum compaction factor. As it has been 
shown in Ref. the function Y satishes the equation 

dY 


uY + iQs 


00 


^Qav[p{0)Y{0,u) 


F{0, u)Y (6*, u)du 


' —OO 

/*oo 


= 2 


q{0'-0)eyip[i{C,-iy){0-0')\d0' / F(A')y(A', 0') dw' 


( 6 ) 


where Qs is synchrotron tune, /AQav is betatron space charge tune shift averaged on all 
variables, and q{0) is the reduced wake function which is connected with usual wake held 
function Wi(z) by the relation: 

roRNWi{-R0) 


q{0) = 


Svr/^yQo 


(7) 
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with ro = e^/mc^ as classic radius of the particle, R as the accelerator radius, N as the bunch 
population, (3 and 7 as normalized velocity and energy. Dehnition of the wake function 
and numerous examples can be found in Ref. [^. More often than not, this function is 
negative (resistive wall wake can be mentioned as the best known example). Just such a 
case was considered in the most of published papers. However, positive wakes are known as 
well, for example the held created by heavy positive ions in a proton beam 1^. It will be 


shown in Sec.VII that the wake sign is especially important when space charge is included 
in the consideration. Alternating wakes are possible as well (resonator models [^) but this 
case is beyond the scope of the paper. 

It is assumed in any case that q{6) is rather short range function so that the wake held 
cannot reach subsequent bunches or turns. Constant wake q = qo will be used to start with, 
and more general cases will be considered in Sec. 6 . Besides of these, space charge is not 
taken into account in 3 nearest sections. With this reservations, Eq. ( 6 ) obtains the simpler 
form 

/■°° 

uY + iQs—= 2q(iexp{iC^9) exp{-iC^9')d9' F{A')Y{A', (p') du' ( 8 ) 

d(p Je J-oc 

where — u. Note that, in most cases, the addition u is rather small to identify (i, 

with normalized chromaticity given by Eq. (5). 

III. EXACT SOLUTIONS FOR HOLLOW BUNCHES 


The hollow bunch model is characterized by the distribution functions 


( 9 ) 


where Aq = ^0 is synchrotron amplitude of any particle and, simultaneously, the bunch 
half-length. The model was repeatedly investigated among others in frames of separated 
or truncated multipole approximations (see e.g. [13|). However, its special interest is just 
that the exact solution can be obtained without any similar assumptions. Therefore corre¬ 
sponding results will be substantially used in this paper to develop adequate approximate 
methods with more realistic distributions, and to control their accuracy. 

The only amplitude A = Aq is essential in this case. Therefore Eq. ( 8 ) with = 0 can 
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be reduced to the one-dimension equation for new function Y^cf)) = Y (Aq, 0 ): 

z/F(0) + z QsY\ct)) = ^ r(00 d0' (10) 

where [0] is the periodic polygonal function of period 27r taking the value [0] = |0| at 
|0| < TT. It is convenient to separate even and odd parts of the function presenting it in the 
form F(0) = 1O_(0) -I- F_(0) with 1O_(0) = Y'+(—0) and l^-(0) = —Y_{—(f)). It is easy to 
see that the even part satishes the equation 

QlYl(4,) + ^ r Y44,')d,l,' (11) 

^ Jo 

Restricting the consideration to the region 0 < 0 < tt, one can transform it to ordinary third 
order differential equation 

Y:'W) + ^Yi(<p) = ^YA<P) ( 12 ) 

his ^his 

with the boundary conditions 


n(0)=0, QlY1(0) + v'^Y\{0)=a, y:(>r) = 0. 


(13) 


Eq. (12)-(13) have been obtained hrst in Ref. [131] where graphical solution id displayed as 
well. The analytical solution is represented below in the standard form: 




(14) 


where Ai _3 are the roots of the algebraic equation: 


\3 I \ _ _ r\ 

X + X — g, g — -(15) 

The constants have to be determined trough the boundary conditions (13). The substi¬ 
tution provides a series of linear uniform equation for Ci _3 which is solvable if corresponding 
determinant is equal to zero. It results in the dispersion equation for the bunch eigentunes v. 


Ai(A 2 — A3)(1 — A2A3) exp( 7 rz/Ai/Qs) -|- 

A2(A3 — Ai)(l — A3A1) exp(7rz/A2/Qs) + ( 16 ) 

A3(Ai - A2)(1 - AiA2)exp(7rz/A3/(55) = 0. 

Thus, the following steps are evident for handling the problem: (i) to hnd the roots Ai _3 
of the hrst Eq. (15) with arbitrary chosen parameter g-^ (ii) to use them for the solution 
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FIG. 1: Real eigentunes of hollow bunch without chromaticity (exact solutions). There are complex 
roots at \qo\/Qs > 0.567 (not shown). 

of Eq. (16); (iii) to substitute the found roots v in second Eq. (15) obtaining dependence 
by exclusion of the parameter g. 

Generally, this way is applicable for all roots including complex ones. However, in this 
section we will restrict our consideration to real roots only, postponing discussing of complex 
eigentunes to next parts where they will be analyzed with chromaticity taken into account. 

Eq. (16) has inhnite number of solutions some of which are plotted in Fig. 1. It is seen 
that V ~ mQs at small go which case has to be treated as independent oscillations of 
different multipoles. Because these tunes are real numbers, instability is impossible at small 
wake and zero chromaticity. Complex roots appear for the hrst time at |go| > 0.567as 
a result of the coalescing of multipoles m = 0 and m = ±1, dependent on sign of the wake. 
Presented inequality should be treated as the threshold of the lowest TMCI mode. Higher 
TMCI modes are possible too being caused by a coalescence of higher multipoles. They have 
essentially larger thresholds as it is shown in Table I (only positive wakes are considered in 
the Table because the picture is symmetric). No other coalescences and TMCI appearances 
have been observed at the calculations. 

Hence the modes |m| > 1 have a minor interest from point of view of the instability 
threshold. Therefore our goal should be an investigation of the lower modes presented by 
red line in Fig. 1, using realistic distributions and wake functions, with chromaticity and 
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TABLE 1: Threshold of different TMCI modes 


Coalesced multipoles 

(0-1) 

(2-3) 

(4-5) 

(6-7) 

(8-9) 

(10-11) 

(^o)t/iresh/ Qs 

0.5671 

3.459 

7.366 

11.894 

16.871 

22.198 


space charge taken into account. 


IV. INTEGRAL EQUATIONS FOR MULTIPOLES 


General solution of Eq. (8) can be presented as the Fourier series: 

V{A, 0) = ^ Ym{A) exp(im0) (17) 

m 

which is just expansion over the multipoles. Its substitution to Eq. (8) allows to get set of 
coupled integral equations for the functions Ym{A) [^: 

poo 

(n - mQs) Y^A) = 2nqo V / K^,n{A A')Yr,iA') F{A')A' dA' (18) 

n Jo 

with the kernels 

^,.Tr,ie'/A')d9' 


2 r . T^{eiA)de 


-A' 


Km,n{A, A') = — / exp«^6') 


TT" J-A 

where Tm{x) = cos(m arccosx) are Chebyshev polynomials, 
6a' = 9 at A^ > 0^, and 


,_ , exp(—i(^9' , _ 

Je^, VA'2 - 


(19) 


9A’ = A' X sign(6') at A'^ < 9'^. 


The assumption = 0 is used to start with but chromaticity will be added hereafter. 

As it was established in previous section, the beam instability is governed mostly by 
the multipoles m = 0, ±1, whereas other ones are of a minor importance. Corresponding 
kernels are: 

A"m,m(A, A') = iLm,-m(A, A') = (20) 


-^o,±i(A, A') — —iL±i^o(A', A) — 


^ \t A'Z LA A 


TtM' V A2 -02 


( 21 ) 


(note that Eq. (20) is valid with any m at C,^ = A). Therefore equations for three mentioned 
multipoles can be written in the form 


vYa(A) = 2nq„ /“{ Y„{A') + K„AA, A) [y_,(.4') + Y,{A')\]PiA)A dA 


( 22 ) 
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(i^ T Q,)K±i(/l) = -2ir* / Koj(A', A)Yo(A')F(A')A'dA' 


(23) 


Excluding Y±i one can obtain the ordinary integral equation for the function Yq{A)\ 


uYoiA) = 2nqo / YoiA')F{A')A' dA' - / K(24, A')YoiA')F{A')A' dA' (24) 

with the kernel 

poo 

K(A4l')= / Ko,i{A,X)Ko,M',X)F{X)XdX (25) 

Jo 

Multiplying Eq. (24) by Yq{A)F{A)A and integrating over A, one can get the simple equa¬ 
tion for the eigentunes v 

u - n / n 2 \ 

(26) 


u — q 


u -= -a go 

u 


with the parameters which do not depend on u: 


q = 


2ngolJo°°yi>{A)F(A)AdA]- 
j;^Y'MF(.A)AdA 


and 


a2 = 


47r2 K(24, A') Yo{A)Yo{A') F{A)F{A') AA' dAdA' 


(27) 


(28) 


ri7(A)F(A)AdA 

These parameters are not too much sensitive to choice of the function Yo(A), so that some 
approximate solution of Eq. (24) can be used for the estimation. The function Yq = 1 
is a simple and rather reasonable option because (i) it satishes Eq. (24) without coupling 
presenting the lower radial mode at go 1, (ii) and it is one of the exact solutions of Eq. (12) 
for hollow bunch with any gg. True, higher radial modes are excluded from consideration by 
this choice. However, it does not matter in the case because it is known that these modes 

a.e.uch.„.esU6.eSa. 


With this choice, the relation g = go follows from Eq. (27), and Eq. (28) also obtains a 
compact form (see Appendix): 


a 


= — F{A)A^ dA / p(A cos 0) sim 0 


(29) 


Jo '-Jo 

Kesnlts of calculated by this formula are collected in Table II for several distributions. 
Somewhat surprising fact is very slight dependence of the parameters on the bunch shape, 
even for so far apart models as the hollow and the Gaussian bunches. It is conhrmed by 
Fig. 2 where solutions of Eq. (26) are plotted against the wake strength, and the images of 
different distributions are also indistinguishable. The data taken from Section III are also 
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FIG. 2: Eigentunes of hollow bunch without chromaticity (3-modes approximation). Red and blue 
curves present real and imaginary parts of the tune. Solid lines are obtained by solution of Eq. (25), 
and some data are transfered from Eig. 1 being presented by circles. 


transferred in this graph being presented by the dark circles, and demonstrating absolnte 
agreement of the resnlts obtained by so different methods. 

The eigentnnes shown in Fig. 2 are real nnmbers at low wake. However, they obtain an 
imaginary addition (bine lines) at rather large qo, which is the TMCI threshold in the case. 
Corresponding values are presented in the last line of Table II being about 0.57 in all the 
examples. In agreement with definition given by Eq. (7), it allows to represent the TMCI 
threshold in usual terms 


roRNlWil , , 

, ' ^ ' > 0.577r ~ 1.8 30 

which expression almost does not depend the bunch shape. The result coincides very well 


with the known expressions 


flK, 


li. 


TABLE II: Parameters a and TMCI threshold of different bunch shapes 



Hollow* 

Hollow 

Rectang 

Parabol 

Gauss 

a 

N/A 

0.4053 

0.4082 

0.4075 

0.4020 

11 threshold /Qs 

0.5671 

0.5689 

0.5672 

0.5676 

0.5708 
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V. THREE-MODE MODEL WITH CHROMATICITY 


As it has been shown in previous section, the function Yq = 1 is rather good approxi¬ 
mation to describe the single bunch instability near threshold at zero chromaticity. Besides, 
the relation 

i*7r 

l±i(A) oc A / p(A cos0) sin^ 0 (31) 

Jo 

follows from Eq. (21) and Eq. (23) in this approximation. Important thing is that the 
integral in this expression moderately depends on the amplitude at realistic distributions. 
For example, it does not depend at all for the rectangular bunch, and has a variation not 
more of 25% for the parabolic one. Therefore the approximation Y±i oc A looks rather 
reasonably for this case. It means that all above presented results could be obtained using 
the pattern solution 

Y{A,(P) = l + Ce9 + CuU (32) 


with indehnite constants Cg and Cu- Conhrmations of this statement will be furnished later. 
However, the main thing is that this model paves the way to extend the theory by including 
chromaticity, space charge, etc. Chromaticity is the hst point which will be applied in this 
section. 

Substitution of Eq. (32) to Eq. (8) results in 

poo 

u + {iyCg+iQsCu)0 + {iyCu-iQsCe)u = 2qoexp{iC6) / p{6’){l + Ce9') exp{-i(6') d9' (33) 

Je 

The relation Cu = iCgQsjv follows from this immediately. Two more equations can be 
obtained by multiplication of Eq. (33) by p{9) or p{9)9 with subsequent integration over 9. 
Then, excluding parameter Cg, one can get required dispersion equation for the eigentunes 
V. It is represented below for the case |C^o| <~ 0.5 which assumption allows to estimate 
effect of chromaticity without excessively bulky expressions: 

v-qQ{l-iax)( \ 2 /^ 

-2-( 1 ^( a-^ 1 (34) 

where y = 2\/2(,yag, and following designations are applied: 


a 


2 

e 



p{9) 9^d9, 


(35a) 



c^e 


p{9)9d9 / p{9')d9', 


(35b) 
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1 


1 


9 


(35c) 


/9 = 


a: 


\V2J-1 


p{e)ede / p{e')e'^de' 


Because ag is rms bunch length, x can be treated as betatron phase advance caused by 
chromaticity in the entire bunch (it is really true for the hollow bunch of length 29q when 
ce = 9q/\/2). Other parameters are represented in Table III for several distributions. Com¬ 
parison with Table II let us to conclude that Eq. (26) and (34) coincide not only formally 
but also actually at y = 0, because the difference of coefficients a is negligible. It can be 
concluded as well that the dependence of the eigentunes on the bunch shape is very weak 
when the chromaticity is also included. The statement is conhrmed by Fig. 3 where com¬ 
plex solutions of Eq. (34) are plotted against the wake strength at different chromaticity, 


TABLE III: Dispersion equation parameters of different bunches 



Hollow 

Boxcar 

Parabolic 

Gaussian 


9112 

0 o 73 

el/f> 

Any 

a 

0.405 

0.408 

0.407 

0.400 

/3 

0.135 

0.123 

0.113 

0.100 
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<: 







FIG. 3; Eigentunes with chromaticity: left - hollow bunch, right - Gaussian distribution. All the 
curves are odd functions of the wake strength. Real parts do not depend on sign of chromaticity, 
imaginary parts change sign. Some curves are indistinguishable because the modes starting from 
the point u = —Qs almost do not depend on chromaticity. 
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for hollow and Gaussian bunches. Note that only positive wakes and chromaticities are 
presented in the graphs because all curves have following symmetry properties: (i) they are 
odd functions of qo/Qs', (h) real parts of the tunes do not depend on sign of chromaticity; 
(iii) the imaginary additions reflect specularly with respect to the line z/ = 0, when the 
chromaticity change sign. It is seen that the instability has no threshold with chromaticity, 
and TMCl appears against the head-tail background without a pronounced demarcation 
line. In particular, the head-tail and TMCl contributions are comparable at |go| — ^-7Qs 
and l^oCI — 0-5- It should be noted in addition that no sign of chromaticity can prevent 
instability of all bunch modes. 

With accuracy of several percents, solutions of Eq. (34) at go < Qs can be presented in 
the form 


u ~ 


go -h mQs - igox(« 



mQs - igox(« + /^) 




a — 




(36) 


where m = ±1. In particular, it provides correct TMCl threshold without chromaticity, 
and leads to well known formulae for the head-tail modes at go -C Qs- hi the last case the 
solutions for hollow bunch can be reduced to the form: 




niQs -h go5m,o + 


^ 2 ( 4^2 _ 


(37) 


which expression is valid with any m |7|]. Analyzing Table III, one can add that this result 
almost does not depend on the bunch shape, at least for lowest radial modes and multipoles. 


VI. REALISTIC BUNCH WITH ARBITRARY WAKE 

It would be beyond reasons to treat rectangular wake as an exclusive case. On the 
contrary, Eq. (32) can be applied as an approximate solution of general Eq. (6) to look for 
the eigentunes of a bunch with arbitrary wake function. Subsequent transformations are 
described just after Eq. (32) and result in the dispersion equation like Eq. (34) or (26): 

2a2ge/^ = -aWef (38) 

with the coefficients 

noo POO 

qef = 2 / q{e)de / p{e' - e/2) p{e’ + e/2) de' (39a) 

Jo J —oo 
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«! = 


o’eqefV^ Jo 


q{e)edex / p{9' -e/2)p{e' + e/2)d9' 


(39b) 


tt2 — —n - / 

mef Jo 


q{9) d9 / p{9' - 9/2) p{9' + 9/2) (r - ^74) d9' 


(39c) 


where oq is rms bunch length given by Eq. (35a), and q{9) = q{9) exp{—i(i,9). Approximate 
solutions of Eq. (38) can be presented in the form like Eq. (36): 


ge/(l + a2) +mQs . 
u ~ — -± 


ge/(l - tt 2 ) - mQs 




(40) 


with m = ±1. Although parameter y does not appear in the expression, chromaticity is 
still presented here being included in the functions q{9) and gef- However, next consideration 
will be restricted by the case of zero chromaticity: = 0, q{9) = q{9). Then 0 : 1^2 are real 

numbers, and the instability can appear only in the TMCI form with the threshold: 

Qs Oi2 




thresh 


a. = OLi — 


l + 2a’ 2 

For Gaussian bunch with dispersion cr^, used parameters obtain the forms 


(41) 


= 


«! = 


02 = 


Jo 

1 

2^/^ajq, 


exp 


9^ 




q{9) d9 


exp 


ef JO 

1 


4cr2 


qi9) 9d9 


9^ 

2al 


exp 


02 


Aal 


q{9) d9 


(42a) 

(42b) 

(42c) 


A^/7iaeqef Jo 

Note that q^f = qo, ai = a, 0:2 = 0 at constant wake q = go- Two more examples are 
considered below. 


A. Short rectangular wake 


Gaussian bunch with a rectangular wake of restricted length 0^ is considered in this 
subsection as first example (similar wake can be created e.g. by a pair of strip-line BPM 
|l^). Eq. (42) gives in this case 


ge/= goerf(a:), «! = 


1 — exp(—a:^) 
\/^erf (x) 


0(2 = 


a;exp(—a;2) 
27 ^erf (x) 


(43) 


14 




















X = 0w/(2a) 


FIG. 4: Parameters of Gaussian bunch with short rectangular wake. Brown line describes an 
effective weakening of the wake in a short bunch, red line represents TMCI threshold. 

where x = 9^/{2a0). These functions are plotted in Fig. 4. Threshold value of go is shown 
as well being determined with help of the expression 

\Qo\thresh _ 1 (AA) 

Qs [1 + 2a(a;)] erf (x) 

It is seen that a shortening of the wake results in a rise of the threshold which becomes 
especially noticeable at <~ 2(T . 

B. Resistive wake 


Resistive wall impedance is the most general and important source of transverse insta¬ 
bilities in circular accelerators. Corresponding normalized wake function is 

tqR^N 


t n\ 

= To' 


Qrw — 


(45) 


27r/3'yQob^ V Rac 

where b is the beam pipe radius, and Uc is the pipe wall conductivity (see e.g. Q). With 
this wake, integrals in Eq. (42) are representable in terms of gamma functions: 


Qef = 


gr^r(l/4) 1.4464 g^ 


Oil — 


r(3/4) 


= 0.2390, 


0^2 — 


(46) 


s/^ \/2r(l/4) 

Threshold value of the effective and usual wakes can be found then with help of Eq. (41): 
\Qef\thresh = 0.739 Qg, \qrw\thresh = 0.511 Qs\/^g- Therefore the resistive wall TMCI thresh- 
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old in usual terms is 


TqR -/Vtferesfe / ^ _ Q Cl ('47') 

27rP^QoQsb^y a,a,~ ' ^ ’ 

where standard rms bunch length cr^ = aeR is used. However, it is necessary to take into 

account that Eq. (45) is valid only at R6 >~ 6/7 when the wake reaches a maximum. 

Therefore sufficient condition of applicability of Eq. (47) is 3> b/'j. 

Another restriction comes from the fact that the resistive wake has a long and slowly 
decaying tail. Therefore it can impact not only next bunches but also itself by the succeeding 
turns. These multibunch/multiturn collective effects should be included in a comprehensive 
investigation of resistive wall instability. However, this point is beyond the scope of the paper 
where only single bunch effects are examined. Nevertheless it can be noted that presented 
results give a possibility to estimate a relative danger of the effects by a comparison of the 
contributed tune shifts. Indeed, TMCI of a single bunch appears at |i^|tmci — 0.77 Qg as it 
follows from Fig. 1-3, and corresponding bunch population is determined by Eq. (47). The 
collective modes i.e. long term tune shift with this intensity is 121 


klLONG — 0.51 Qs 


(, _ (2|fc-(3o|77 

27 rR\k - Qol V ) 


(48) 


where h is number of bunches, and k is the collective mode number (the mode can be 
unstable at A; > Qq). Taking \k — Qq\ = 0.25 and (3 = 1 as a typical example, we can 
compare these long term and TMCI effects as the ratio of corresponding tune shifts: 

0.35' 


^^LONG 


l^TMCI 


2 h- 
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2ttR 


(49) 


With great probability this value is < 1 or even -Cl at h = 1, that is the multiturn effect 
of a single bunch is typically small or even negligible in comparison with TMCI or head-tail 
instability. However, the collective effects are more dangerous in a multibunch machine with 
h 3> 1, haz ~ R- Of course, more detailed analysis is required at intermediate cases. 

VII. SPACE CHARGE EFFECTS 


Bunched beam instability with extremely large space charge (AQ S> Qs) was considered 
in works 1^12]. The most remarkable phenomenon is a pronounced asymmetry of the 
curves with respect to the wake sign which effect has been Erst shown in Ref. Hj]. It is 
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FIG. 5; Eigentunes of a rectangular (“boxcar”) bunch with ultimate space charge AQo S> Qs- The 


graph is taken from 


lH; n = m{m + l)/2. 


illustrated by Fig. 5 taken from the quoted article where a rectangular (“boxcar”) bunch 
with constant wake was explored at space charge betatron tune shift AQo ^ Qs- It is seen 
that TMCI appears only with the positive wake satisfying the instability condition 

0.5 Q2 


iSlQthresh — 




(50) 


which drastically differs from the conditions presented by Fig. 1-3, and Eq. (30) of this 
paper. It is needless to say that investigation of the effect at AQo ~ Qs is the only way to 
resolve the problem by a joining of these conflicting pictures. It turns out that the 3-modes 
model described by Eq. (32) bridges these ultimate cases providing a general form of the 
lowers eigenmodes over a wide range of parameters. 

It is easy to verify that substitution of Eq. (32) to general Eq. (6) results in 


v -\- {i'Ce-\-iQsCu) ^ T { [ {^Qavp{(^) + Cu — iQsCg ]} u 

poo 

= 2qoexp{i(6) / p{6'){l + CgO') exp{—i(6') dO' 


(51) 


The only distinction of this expression from Eq. (33) without space charge is the addition 
proportional to AQavp{0). We consider in this section a rectangular bunch which density 
does not depend on longitudinal coordinate so that QavP is the constant value AQq which 
coincides with incoherent space charge tune shift averaged over all transverse coordinates 
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12| . Therefore relation between coefficients of Eq. (51) obtains the form Cu = iCgQs/{v + 


AQo) instead of Cu = iCgQs/v- As a result, all subsequent relations hold true with the 
replacement of Ql/i' on Ql/{v + AQq). In particular, Eq. (34) with constant wake function 
and chromaticity obtains the form 

z/- go(l - *ax) Ql 


V — 




^X 


V + AQo 

In the “head-tail” limit, that is at |q'o| Qs, approximate solutions of the equation are 

AQq AQo 


(52) 


z/Q = go(l-*ax), ^±i-M Ql + 


+ 


(53) 


Thus space charge does not affect zero mode at all and does not change growth rate of the 
modes m = ±1 in this “head-tail” approximation. 

Another situation arises at Qq ~ Qg or go > Qs where TMCI can arise. This case is 
illustrated by Fig. 6 where the eigentunes are plotted against the wake strength at different 
space charge, but without chromaticity. Effect of the wake sign is seen very clearly in this 
graph; space charge propels the TMCI threshold to the centerline at go > 0, and away 
from it at go < 0. Corresponding dependence is shown quantitatively in Fig. 7 where the 
thresholds are presented separately for positive and negative wakes and supplemented by 
appropriate analytical formulae. There is very good agreement of these results with Fig. 5. 
For example, TMCI threshold of positive wake is 


(?o)z 


0.57QI 


-)■ 


0.57 Q 


(54) 


Qs + AQo AQo 
what is very close to the estimation given by Fig. 5 and Eq. (50). 

However, space charge tune shift raises the TMCI threshold of negative wakes. For 
example, Eq. (47) for resistive wall TMCI threshold obtains the form: 

TqR Nfhresh 


27r/3'yQoQsb^ 
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-0.51 1 + 
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Qs 


+ 




(55) 


Joint effect of space charge and chromaticity is illustrated by Fig. 8 where dependence 
of the eigentunes of rectangular bunch on the wake strength is presented at AQo = Qs and 
various chromaticity. 
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FIG. 6: Eigentunes of rectangular bunch with space charge (chromaticity is turned off). Used 
ratios AQq/Qs are: 0 (maroon), 1 (red), 2 (green), and 3 (blue). Solid/dashed lines represent 
real/imaginary parts of the tune. 



FIG. 7: TMGI threshold of the rectangular bunch against the space charge tune shift. Positive 
and negative wakes are presented separately, no chromaticity. 
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FIG. 8: Eigentunes of rectangular bunch with AQq = Qg and chromaticity. Solid/dashed lines 
represent real/imaginary parts of the tune. Some red and blue solid lines coalesce because real 
part of the tune is even function of chromaticity. 


A. Discussion 

The most unstable TMCI mode with space charge was considered in Ref. {q] where ex¬ 
pansion of eigenfunctions in terms of azimuthal and radial modes has been applied. Fig. 1 
of the paper gives an example of negative wake which is twice the size needed to produce 
instability without space charge: kF/|kFthresh(0)| = —2. The results depend on number 
of the basis modes which is characterized by the number m^ax- According them, the in¬ 
stability disappears at rather large AQsc/Qs, and threshold values of this parameters are: 
0.85 at TUmax = 1 (3 multipoles), 0.5 at rumax = 5 or 10 (21 multipoles and up to 5 radial 
modes in the last case). My Eq. (52) provide reasonably close parameters of the threshold: 
AQq/Qs = 1 at qo/Qs = —1.14 (point X = 1, E = —1.14 in Fig. 7). 

However, there is a profound disagreement at the larger space charge tune shift. A 
monotonous behavior of the threshold follows from my paper, and this statement does not 
contradict Fig. 1 with m m a x = 1 and 5. However, new region of instability at AQ > 2.2 is 
predicted by the figure with mmax = 10. This problem is pointed but not explained in [9j 
(end of Sec.H). 
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Analytical and numerical investigation of the model with hollow bunch in square potential 
well is also performed in Ref. {qI- Being compared with Fig. 1, the results agree with option 
^max = 1 better than with multi-modes approach. At least, Fig. 13 obtained by numerical 
solution of differential equation demonstrates a monotonous behavior of the threshold, with¬ 
out additional instability regions at higher AQq. However, analytical expression for AQx in 
page 10 gives a way of head-tail instability even without chromaticity and space charge. I 
guess this statement is in a conflict with commonly accepted point of view [?, 8, 

Eq. (53) has another appearance and does not suffer from this shortcoming. 


IM- My 


My results correlate well with the multi-particle simulations presented in Ref. m . For 
example. Fig. 1 of this paper resembles my Fig. 6 and allows to hnd TMCI thresholds of 
square bunch at AQq/Qs = 4. According it, the relative wake strength is: IFthresh/lfhol — 
for negative wake, and 0.2 - 0.3 for positive one where Wq is the threshold value for no 
space charge (it is difficult to get more exact numbers from the plot). In terms of my paper 
hFthresh/|hFo| = (<?o)thresh/(0.57Qs) wffich value should be —8.6 or 0.2, according to my Fig. 7. 
Another example is provided by Fig. 3 [l^ where threshold value of the wake is presented 
as a function of AQ/Qg. The curves for square bunch coincide with my Fig. 7 not only 
in shape but also quantitatively. Indeed, it follows from Fig. 3 that IFthresh/|lTo| — —6 
at AQo/Qs = 3. In terms of my paper, it means qthresh/Qs = —6 x 0.57 ~ —3.4 while 
the value -3.3 follows from Fig. 7. Positive wake thresholds are in a good consent as well 
being presented in these hgures. However, the agreement is not so close for smooth bunch. 
According to Fig. 3 16|, thresholds of the smooth and square bunches have a similar behavior 
at AQsc/Qs < 3. The similarity is especially obvious if the averaged across the bunch value 
is used as the argument (AQo = 0.43AQsc for this distribution). However, the results come 
apart at larger tune shift because non-monotonous behavior of the smooth bunch is shown 
in Ref. jl^ . 


VIII. CONCLUSION 

The theory of a single bunch transverse instability is advanced in the paper by devel¬ 
opment of 3-modes model for the most unstable bunch modes. The dispersion equation is 
presented in form of 3rd order algebraic equation which includes chromaticity and space 
charge, and can be used with any bunch shape and wake held form. 
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The known TMCI and head-tail instability appear in the theory as the limiting cases. It is 
shown that a distinct boundary between them exists only at zero chromaticity representing 
the TMCI threshold in the case. Generally, the TMCI appears more or less smoothly against 
the head-tail background without a pronounced demarcation line. 

The results depend very slightly on the bunch shape so that simple bunch models can 
be successfully used to analyze the stability limits. For example, difference of the TMCI 
thresholds is less of 1% for so far models as hollow and Gaussian bunches, if they have the 
same rms length and space charge is negligible. 

In contrast with this, the tunes essentially depend on the wake form. Several cases are in¬ 
vestigated in the paper including arbitrary rectangular and resistive wall wakes. Comparison 
of the single bunch and multibunch/multiturn effects is realized in the last case. 

Space charge tune shift is included in the consideration at arbitrary relation of the shift 
to the synchrotron tune. It is shown that the space charge effect depends on the wake sign: 
it increases the instability threshold if the wake is negative, and decreases it at positive 
wakes. Simple analytical formulae are presented for the instability threshold and growth 
rate. They coincide well with the known expressions in the limiting cases though generally 
there are some divergences which are discussed in the paper. 
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IX. APPENDIX: DERIVATION OF EQ. (29) 


It follows from Eq. (25) and (28) at Iq = 1, 


a2 = 87r=^ / F{A) AdA 


Ko,i{A',A)F{A')A' dA' 


uo 


The substitution of Kq i from Eq. (21) results in 
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where 6a = minj^, A}. Changing sequence of the integrals obtain 
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The last integral is p{d)/2 so the expression is reducible in the form 
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which can be transformed in Eq. (29) by the substitution ^ = cos </>. 
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